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ABSTRACT 

For use in calculating higher-order coherent- and squeezed- state quan- 
tities, we derive generalized generating functions for the Hermite polyno- 
mials. They are given by J2™=o z ^ n+k Hj n+ k{x)/(jn + k)\, for arbitrary 
integers j > 1 and k > 0. Along the way, the sums with the Hermite 
polynomials replaced by unity are also obtained. We also evaluate the ac- 
tion of the operators exp[a^ (d/dxY] on well-behaved functions and apply 
them to obtain other sums. 
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1 Introduction 



In calculations involving coherent and squeezed states, infinite sums involving Hermite 
polynomials are extremely useful. Even so, knowledge of such sums is limited. For 
example, if we define the generalized generating functions for the Hermite polynomials 

as 

oo ^jn+k J_ 



G{j, k) = J2 U " ,+I - U) 



for arbitrary integers j > 1 and k > 0, it is only the specific generating functions 



G(1,0) = exp[-z 2 + 2xz\ , 
G(2,0) = exp[-z 2 ] cosh(2x;z) 
G(2, 1) = exp[-z 2 ]smh(2xz) 



(2) 
(3) 
(4) 



that are well-known [T]. Physically, Eqs. (|2])-(^|) are very important as they describe 
the decomposition of ordinary, even, and odd coherent states, respectively, into the 
harmonic-oscillator number states ||. 

Motivated by this situation, in this paper we will show that the general generating 
functions |J G(j,k) are given by 



1 * exp [-z 2 e i4nl/j ] exp [2xze i2nl/j ] 

111 [ e i27rZ/j]fc 



Along the way we will also obtain the sums 



S(j, fc) = E 



00 gjn+k 



1 j 

E 



exp[2e 



i2wl/j] 



f (jn + k)\ 3 ti [e i2 ^} k 
We will use the same techniques to evaluate 



L = exp 



dx 



g(x) 



(5) 



(6) 



(7) 



and apply it to obtain further sums involving Hermite polynomials. 
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2 The Sums S(j, k) 

First note the special cases 

5(1,0) 
5(2,0) 



cosh z = - 
2 



e z + e 



(8) 
(9) 



A few other sums are known. In particular, 5(3, 0) can be obtained from the known 
formula for 5(3, 1): § 



3 M = [x) 3 M=[£ 



I e *_2 e -^ 2 cos U^ z _ + 1 
3 3 V 2 3 



(10) 



With a little algebra, the result can be written in the form 



1 



5(3, 0) = - (exp[ze i27r/3 ] + exp[ze i4n/3 ] + exp[z] 



11^ 



But from Eqs. 



and nj) we can now guess what the answer for S(j, 0) is. 



It is 



1 



ze 



iZirl/j] 



J 



(12) 



i=i 



The phase factors in the exponentials are the jth roots of unity. Therefore, in the 
sum of the Taylor series expansions of the exponentials, after the z° term, the terms 
proportional to z, z 2 , . . ., z^~ x cancel to zero. The next non-zero term is the z^ term. 
Similarly for the higher-order terms. 

This argument can be made precise by using the geometric series 



i=i 



r — 1 



(13) 



Taking r = exp[i2ir/j], r = exp[i47r/j], r = exp[i67r/j], . . ., gives the result. 
Then, S(j, k) of Eq. @ is obtained by observing that 

k 



S(j,k)=(Jdz^ S(j,0). 



(14) 
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3 The Generating Functions G(j, k) 

Now we are ready to calculate the generating functions G(j, k). In Eq. ([!]) for G(j, k), 
substitute the definition flll 



Jn 

H n (x) = (-l) n exp[x 2 } — exp[-x 2 ] 



One then has 

G(j, k) = exp[x 2 ] 
which from Eq. (0) is equal to 



oo / y d \ jn+k 

-y^ v f dx> 
h Un + k)\ 



exp| — x 



(15) 



(16) 



k) = exp[x 2 



1 ^ exp[-ze i2wl/jA 



dx J 



exp[— x z 



(17) 



But now using the relation (see below) 

d 



cxp 



dx 



g(x) = g(x + a) 



one obtains the final result for G(j, k) given in Eq. (|5|). 



4 The Operators exp[a J \d / dx 



The techniques we have used above can be applied to other situations. For example 
we may ask if the result of Eq. ([18]) , 

d 



I\ = exp 
can be generalized to j > 1: 
Ij = exp 



dx 



9{x) = 2^ = 9{x + a) , 



n=0 



111 



dx 



, ~ a jn g^ n Mx) 
9{x) = }^ 



n=0 



rv. 



First consider the special case j = 2: 



~ a 2n g^ n \x) 
h - 2^ 



n=0 



111 



(19) 



(20) 



(21) 
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Using the doubling formula 



r («+l/2) = -' /! !gj, (22) 



one obtains 

00 {Q„\2n 



= ^|Sf r( " + 1/2)9M(x) 

2 ^£^r(!l±i)[ 9 H ( , ) + ( _ ir9 H (:c) ] 



ju b n=0 lb. 



2^/2 y s 1 / 



dS i~[ff(a: + 2a* 1 / 2 ) + g(x - 2as 1 / 2 )} . (23) 



2tt 1 /2 Jo sV2 

In the above, the use of (—1)™ allowed us to cancel out the odd derivatives as we 
changed the sum from second-order terms to first-order terms. The use of the integral 
representation of the T function then led to the final result. 
Note that by the change of variables 

s = (y- xf/Aa 2 , (24) 

our result can be put in the form of the known result || 

1 ' w '"'" 9(y)dy- (25) 



POO 


{y-x) r 


/ ex p 

/ — 00 


4a 2 



z (4a%)V2 

(This formula will be useful in obtaining further sums in the next section.) 

The way to generalize to arbitrary j > 2 is now, in principle, clear. First one uses 
the multiplication formula (which is the generalization of the doubling formula) in 
the form 



J 



+1/2 Ul' = \T(n + k/j) 



n\ (2vr)0'- 1 )/2 {j n )\ ' ^ 

Then, just as before, the use of exp[i2irl / j] will allow the cancelation of undesired 
terms in the complete sum over n. Here those terms are the g^ n ^. When all is done, 



the end result is 

1 r°° dxi e~ xi f°° dx 2 e~ X2 r°° dxj-i e~ Xj ~ 1 

3 = (27r)«-i)/2ji/2 J x i/j J ~JTi •••./„ JfWi 



x 2 " JU a;j_i 



J29 (x + ja(x lX2 . ..x^) 1 '^ 2 ^) . (27) 
l=i 

This result, although in closed form, is usually too complicated to yield results 
in terms of elementary functions. For instance, even for j = 3, changing variables to 
Xi = yf yields 

Q r oo c og ^ 

Is = T^y^ I dyx jf rfy 2 e"^ (x + 3ay l2/2 e jM / 3 ) . (28) 

But there is also another point that must be mentioned. In Ref. |J, it was 
observed that naive higher-order squeeze operators of the form exp[ z j a ~ z *j a3 ]i 
j > 2, are not Gaussian- integrable operators. Therefore, the convergence of the Ij, 
j > 2, depends on the form of g(x). 

Finally, for j > 3, Baker has observed that using the Hankel contour integral 
for 1/r will lead to a formally simpler answer. Starting from Eq. (|20|), one has 



oo j n i2nln/j [n]( \ 



where [(n/j)]! is not clearly defined for n not a multiple of j. However, because of the 
sum over I, only the n-terms which are multiples of j will sum to non-zero. Therefore, 
[(n/j)]! can be replaced by 1/T(1 + n/j), which in turn can be replaced by the Hankel 
contour integral. Unity in the form of T(n + l)/n! can be inserted, and T(n + 1) 
replaced by its integral representation. Then, a little algebra yields 

C is the contour starting just above the x-axis at (oo,ie), coming in and encircling 
the origin counterclockwise, going out to (oo, — ie), and closing. It is straightforward 
to demonstrate that I\ and I 2 of Eqs. fllTf ) and (|23| ) are special cases. 



5 Discussion 



With our methods, it is possible to obtain further quantities. Sums of the type 

00 z jn+k H- ,Jx) 



n=0 (pn + q)\ 

become amenable to analysis, as well as even more complicated sums. 
For example, one can immediately write 



#(2,0,1,0) = exp [x 2 ] exp 
exp [a; 2 ] 

1 



dx 



exp[— x z 



POO 


(y-x) 2 


/ ex p 

1 —00 


4* 2 



e*v[-y 2 ]dy 



exp 



Az 2 x 2 
1 + Az 2 



[l + 4z 2 ]V2 

This describes the squeezed ground state. Also, one has 

d ' 



#(2,1,1,0) = exp[x 2 



2zx 



[1 + Az 2 f/ 2 

Adding these two yields the known sum E 



dx 



exp 



#(2,0,l,0)exp[-^ 



Az 2 x 2 
1 + Az 2 



~ z n H n {x) _ l + 2zx + Az 2 

hWW [T+4^p exp 



2^2 



Az z x 
1 + Az 2 



(31) 



(32) 



(33) 



(34) 



where in the above [[•]] denotes the greatest integer function. 
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